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1.0  INTRODUCTION 

A  series  of  optical  homodyne  experiments  (Refs.  1  through  4)  during  the  period 
1962-1964  demonstrated  the  feasibility  of  exploiting  the  Doppler  shift  of  scattered  laser 
light  for  measuring  the  velocity  of  moving  matter.  Of  particular  importance  to  the 
aerodynamicist  was  the  demonstration  by  Yeh  and  Cummins  in  1964  (Ref.  4)  that  fluid 
velocities  could  be  inferred  by  measuring  the  Doppler  shift  frequency  of  laser  light 
scattered  by  particulate  matter  entrained  in  the  fluid. 

Since  1964,  worldwide  interest,  intensive  development  efforts,  and  increasingly 
refined  applications  have  resulted  in  significant  improvements  in  laser  Doppler  velocimeter 
(LDV)  hardware  and  techniques  as  well  as  in  a  heightened  awareness  of  remaining 
technological  deficiencies.  In  recent  years,  significant  efforts  have  been  directed  toward 
specifying  and  reducing  the  magnitudes  of  the  systematic  and  random  errors  contributing 
to  the  overall  uncertainty  of  LDV  measurements  (Ref.  5).  One  subtle  systematic  error 
mechanism  not  yet  widely  recognized  is  that  caused  by  not  accounting  for  or  improperly 
accounting  for  axis  nonorthogonality  in  multiaxis  LDV  systems. 

This  error  mechanism  was  recently  encountered  in  a  large  transonic  wind  tunnel 
application  of  a  two-axis  LDV  system  (Ref.  6).  This  system  was  used  for  free-stream, 
pitch-plane  velocity  magnitude,  and  flow-angle  measurements.  Orthogonal  axis,  data 
reduction  equations  were  initially  used;  however,  the  measurement  axes  were 
subsequently  found  to  subtend  an  angle  of  about  89.1  deg.  This  oversight  resulted  in 
velocity  magnitude  errors  as  large  as  0.8  percent  of  measured  value  and  flow-angle 
measurement  errors  as  large  as  0.4  deg.  The  tunnel  test  requirements  demanded  overall 
uncertainties  in  velocity  magnitude  and  flow-angle  measurements  of  less  than  ±1.0 
percent  of  measured  value  and  less  than  ±0.1  deg.  Thus,  the  oversight  errors  were 
obviously  unacceptable  and  were  ultimately  eliminated  by  using  the  equations  developed 
in  this  report. 

The  measurement  axes  of  virtually  all  two-  and  three-axis  LDV  systems  are  either 
inadvertently  or  intentionally  nonorthogonal  to  some  degree.  Also,  the  LDV  axes  may  be 
rotated  with  respect  to  a  reference  set  of  rectangular  Cartesian  coordinates.  Whether  or 
not  special  data  reduction  equations  are  required  to  account  for  axis  nonorthogonality 
depends  upon  the  magnitude  of  the  nonorthogonality  and  the  accuracy  sought  in  the 
particular  LDV  application.  When  special  equations  are  necessary,  it  is  essential  to 
incorporate  into  the  transformation  relationships  the  subtle  point  that  the  LDV  system 
measures  velocity  vector  projections  rather  than  components  when  its  measurement  axes 
are  oblique. 
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At  least  two  previous  analyses  have  dealt  with  nonorthogonal  LDV  measurement 
axes.  The  first  (Ref.  7)  treats  a  special  3-space  case  of  transmitter-receiver  symmetry  and 
is  not,  therefore,  universally  applicable.  The  second  (Ref.  8)  provides  only  the  2-space 
transformation  expressions  relating  two-axis  LDV  measurements  to  reference  axis 
components.  A  general,  comprehensive,  3-space  analysis  of  LDV  measurement  axis 
nonorthogonality  and  rotation  apparently  is  not  yet  available. 

The  purpose  of  this  report  is  to  derive  the  general  3-space  equations  for  velocity 
vector  magnitude  and  spatial  orientation  in  terms  of  oblique  axis  parameters  and  the 
velocity  vector  projections  measured  by  each  LDV  axis.  This  is  accomplished  by  first 
employing  fundamental  tensor  relationships  to  derive  expressions  that  interrelate  vector 
measure  numbers,  components,  and  projections.  These  expressions  are  then  used  along 
with  the  metric  tensor  to  develop  the  desired  equations. 

2.0  VECTOR  MEASURE  NUMBERS,  COMPONENTS,  AND  PROJECTIONS 

Figure  1  depicts  a  three-dimensional  space  simultaneously  covered  by  the  rectangular 
Cartesian  coordinate  system  x1  and  the  general  or  curvilinear  coordinate  system  x*'. 


X3 


Figure  1.  3-space  coordinate  systems  x1,*2,*3  and  xr,x2',x3'. 
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Let  x>  and  x'fbe  functionally  related  so  that  the  transformation 

x1'  =  x’Xx1)*  i'  =  l  2',  3  and  1  =  l,  2,  3  (1) 

is  single  valued  and  reversible.  Let  the  inverse  transformation 

X1  =  x'Cx1'),  i  =  1,  2,  3,  and  i'  =  1  2',  3'  (2) 

also  be  single  valued.  Under  these  conditions  the  variable  point  P  shown  in  Fig.  1  can  be 
uniquely  specified  in  terms  of  either  its  x‘  or  x5'  coordinates. 


Let  r  denote  the  position  vector  of  the  point  P  with  respect  to  0,  the  origin  of  the 
rectangular  Cartesian  coordinate  system.  A  differential  change  in  r  can  be  expressed  by 


dx1 


(3) 


where  the  Einstein  summation  convention  is  used.*  For  unit  displacement  along  the  x1’ 

— i ►  ;  I 

axis,  the  change  in  r  is  tangential  to  x1  and  is  equal  to  3r/3xl  .  The  vectors 

dr 

(4) 


are  the  unitary  vectors  with  respect  to  the  point  P  (Ref.  9)  which,  in  general,  constitute  a 
nonorthonormal  basis  for  all  field  vectors  defined  at  the  point  P. 

The  unitary  set  defines  a  parallelepiped  whose  volume  is  given  by  the  scalar  triple 
product  of  these  vectors  (Ref.  9) 


7?.  —r  ~ r 

U  =  a  |  •  a2'  x  a3 ' 


(5) 


*Summation  over  the  index  for  all  admissible  values  is  to  be  carried  out  whenever  the  same  index  occurs 
both  as  a  superscript  and  as  a  subscript  in  a  single  term.  A  superscript  in  the  denominator  of  a  term  is 
regarded  as  a  subscript  on  the  quotient.  The  present  analysis  is  limited  to  3-space  so  that  the  admissible 
values  are  always,  1,  2,  and  3. 
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A  second,  generally  nonorthonormal,  set  of  basis  vectors  can  be  defined  in  terms  of  the 
unitary  set  and  their  scalar  triple  product  (Ref.  9) 


i'-i  ft.. 


/->  -*■  \ 

“*3  '  1  / 

\a3  '  x  al 7 

-  a  -  c  ( 

(6) 


These  are  the  reciprocal  unitary  vectors.  The  unitary  and  reciprocal  unitary  vectors  are 
seen  to  satisfy  the  condition 


where  5J; 


is  the  Kronecker  delta.* 


(7) 


The  differential  change  in  the  vector  r  can  now  be  expressed  either  in  terms  of  the 
unitary  basis  as 


*T  — ►  ■  ' 

dr  =  aj-dx1 

or  in  terms  of  the  reciprocal  unitary  basis  as 


(8) 


dr  =  aJ  dxj-  (9) 

Analogously,  any  field  vector  V  defined  at  the  point  P  can  be  resolved  into  components 
either  with  respect  to  the  unitary  basis  as 


(10) 


or  with  respect  to  the  reciprocal  unitary  basis  as 

— ►  -it  / 

V  =  aJ  \y  (11) 

-  r 

where  v1  are  the  contravariant  measure  numbers  of  V  and  v are  the  covariant  measure 
numbers.  Equations  (10)  and  (11)  are  statements  of  the  familiar  parallelogram  law  for 
formation  of  the  vector  V  by  vector  addition  of  its  components. 

The  squared  length  of  dr  is  expressed  with  respect  to  the  unitary  set  by 


(ds')2  =  dr  •  dr  =  -  dx 1  dxJ  =  g^-dx1  dx* 


(12) 


S'.  =  1,  i’  =  j' 
S|'  =  o,  i'  j' 
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where  g,'j'  is  a  symmetric  tensor  termed  the  covariant  metric  tensor  (Ref.  10).  This  same 
length  is  expressed  with  respect  to  the  reciprocal  unitary  set  by 


(ds'E 


-  dr 


-> 

dr  = 


ri  dx  -dx  '  =  g1  J  dx  'dx.< 

I  1  &  II 


(13) 


where  g'V  is  a  symmetric  tensor  termed  the  contravariant  metric  tensor.  Tlie  analogous 
expressions  for  a  general  field  vector  V  are 


V2 


£  '  'V 
b,  j 


1  vJ 


(14) 


and 


V2  =  g’  J 


(15) 


A  set  of  unit  vectors  collinear  with  the  unitary  set  is  given  by 


U  .  ■' 

1 


fno  summation  on  D 


(16) 


The  magnitudes  of  the  components  of  V  with  respect  to  this  unit  set  are  given  by 


(V)t,  =  vlv 


v  s,v 


(17) 


The  orthogonal  projections  of  V  onto  the  axes  defined  by  the  unit  set  are  given  by 


-v 

a.  ' 


(Viv 


=\ 


si ' v 


(no  summation  on  1  ') 


V  «j 


V  S;',' 


v 


(18) 


Study  of  Eqs.  (5)  and  (6)  shows  that  the  a,'  and  a1  are  identical  quantities  when  the 
a;'  are  the  orthonormal  basis  vectors  of  a  rectangular  Cartesian  coordinate  system. 
Under  these  conditions 


r .  s  *  —  §  /  /■  —  £ 1  ^  J 

>1  J  yi  j  B  w 


(19) 


and  consequently 

v,*  =  v‘  =  <VV  =  'V,),'  (20) 

However,  v/,  v1  ,  (V),',  and  (Vx)/  are  distinctly  different  quantities  when  the  a,'  are 
nonorthonormal.  It  will  be  shown  subsequently  that  an  LDV  system  measures  the 
quantity  (Vj_y. 


9 


AEDC  TR  79  35 


In  Riemannian  geometric  space  the  element  of  arc  length  ds'  is  given  by  Eqs.  (12 
and  13).  If  there  exists  a  coordinate  transformation 

xr  =  x1  V)  (21) 


such  that 

(ds')2  =  (ds)2  =  dx1  dx^  (22) 

then  the  Riemannian  space  is  said  to  be  Euclidean  (Ref.  11).  In  Euclidean  space,  any 
coordinate  system  for  which  the  elements  of  the  metric  tensor  are  constants  is  termed  a 
Cartesian  coordinate  system  (Ref.  11). 

All  essential  metric  (quantitative  measures  of  length,  angles,  etc.)  properties  of 
Euclidean  space  can  be  completely  specified  by  use  of  the  applicable  metric  tensor  (Ref. 
10).  Specifically,  Euclidean  3-space  covered  by  the  oblique  Cartesian  coordinate  system 
defined  by  the  transmitted  beam  pattern  of  a  multiaxis  LDV  system  is  characterized  by 
a  metric  tensor  with  constant  elements.  Therefore,  it  is  reasonable  to  seek  mathematical 
relationships  for  velocity  vector  magnitude  and  spatial  orientation  in  terms  of  the 
applicable  metric  tensor  and  the  LDV-measured  values  of  the  velocity  vector  projections. 

3.0  LDV  MEASUREMENTS  AND  VELOCITY  VECTOR  PROJECTIONS 


Most  present-day  LDV  transmitters  are  of  the  crossed-beam  configuration  (Ref.  12). 
This  configuration  is  compatible  with  the  reference-beam  receiver  as  well  as  the 
dual-scatter  receiver.  The  expression  for  the  signal  frequency  (Doppler  frequency) 
produced  by  one  axis  of  a  crossed-beam,  reference-beam  LDV  system  is 

fD  =  2^-V.  (t  -tr)  (23) 

where  V  is  the  velocity  vector  of  a  light-scattering  particle  passing  through  the 
crossed-beam  or  probe  volume  region,  kj  is  the  propagation  vector  of  the  illuminating 
laser  beam,  and  ^  is  the  propagation  vector  of  the  scattered,  Doppler-shifted  laser 
radiation.  For  nonrelativistic  velocities,  Eq.  (23)  can  be  put  into  the  form 


2  sin  6/2 


(24) 


where  X  is  the  wavelength  of  the  illuminating  laser  radiation,  6  is  the  angle  subtended  by 
the  two  intersecting  laser  beams,  and  iff  is  a  unit  vector  of  coplanar  with  and  and 
perpendicular  to  the  bisector  of  the  angle  d. 
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The  signal  frequency  produced  by  one  axis  of  a  crossed-beam,  dual-scatter  type  LDV 
is  given  by 

=  17  V  '  (  k  12  -kIl)  (25) 

— ► 

where  k[j  is  the  propagation  vector  of  the  first  of  two  crossed  beams  of  illuminating  laser 
radiation  and  k]2  is  the  propagation  vector  of  the  second  illuminating  beam.  With  9  as 
the  angle  subtended  by  iqj  and  k,2)  Eq.  (25)  can  be  put  into  a  form  identical  to  Eq. 
(24)  which  states  that  an  LDV  system,  whether  reference-beam  or  dual-scatter  type, 
simply  measures  the  orthogonal  projection  of  $  upon  the  u^'  axis. 

4.0  THE  METRIC  TENSOR 

As  previously  stated,  a  general  coordinate  transformation  is  expressed  by 

x'  =  x1  (x1  ) 

The  total  differential  of  x*  is  given  by 


(26) 


(27) 


The  covariant  metric  tensor  is  defined  by  substituting  Eq.  (27)  into  the  expression  for 
scalar  invariance  of  the  infinitesimal  distance  ds  in  rectangular  Cartesian  coordinates  and 
ds'  in  general  curvilinear  coordinates  (Ref.  13): 

Ids')"2  =  gtydx‘  dxJ  =  (ds)2  =  5^  dx1  dxJ  (28) 

From  this  it  can  be  seen  that 


g 


1 V 


(V  dxJ 
dx1'  dx1' 


(29) 


For  the  special  case  of  Cartesian  coordinates,  the  general  transformation  given  by  Eq. 
(26)  simplifies  to 


X 


1 


A 


(30) 


where  A]'  is  a  matrix  of  constants  independent  of  coordinate  position.  Examination  of 
Eqs.  (29)  and  (30)  shows  that  the  elements  of  A]-  must  first  be  determined  before  the 
elements  of  the  metric  tensor  can  be  computed. 
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4.1  COORDINATE  TRANSFORMATIONS  IN  3-SPACE 

In  a  Cartesian  coordinate  system  a  vector  can  be  represented  as  a  directed  line 
segment  connecting  two  points  in  the  space.  Such  a  vector  is  shown  in  Fig.  2  as  the 
directed  line  segment  OP  and  is  identified  as  V.  The  familiar  parallelogram  law  for  vector 


Figure  2.  Components  of  vector  V  with  respect  to  oblique 
Cartesian  coordinate  systems  xr  and  x'\ 


addition  of  the  vector  components  with  respect  to  two  different  oblique  Cartesian 
coordinate  systems  is  illustrated  in  Fig.  2.  Further  study  of  this  Figure  shows  that  the 
orthogonal  projection  of  V  onto  the  x'  axis  is  given  by 

Xl'  +  V2~X2  +br3.X3  =  <Vl'xl  +Cl"2'*2  +Cl'3'x3  (31) 


,  i  t 

where  b,'y»  is  the  cosine  of  the  angle  subtended  by  the  positive  directions  of  the  x1  and 
xi  axes  and  c,"'  is  the  cosine  of  the  angle  subtended  by  the  positive  directions  of  the 
x»  and  x1’  axes.  In  general, 


(32) 


where  b,"j"  and  c,"j'  are  matrix  arrays  of  these  direction  cosines.  Premultiplication  of 
Eq.  (32)  by  b»'v'  yields 


-  b>  ' 


(33) 


where  bJ  ''”  is  the  inverse  of  b,"s”.  A  similar  procedure  produces 


x1  =  bJ  1  c  .  ~  x1 


(34) 
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where  bi  ‘  is  the  inverse  of  bjV  and  c/j"  js  the  transpose  of  Cj"j'. 

For  the  special  case  that  xi  =  xJ,  where  x)  represents  a  rectangular  Cartesian 
coordinate  system,  ”1"  becomes  the  identity  matrix  IJ‘ 

xJ  =  I1’ c  .  ..  Jp  =  AJ,  xJ 

•J  j  (35) 

which  is  in  agreement  with  Eq.  (30).  Equation  (34)  now  becomes 


XJ  =  b’  1  C  -  XJ  =  D5  xJ 
1  J  i 


(36a) 


where 


jVf  []  (  b2'3  ')  2  ]  +  c2'l  (blVb2V  biv)  +  C3l  (bl  Vb2V  b  1  '3  ')  J 


(36b) 


DC  A  '  2  1  (b£'3')2_  *  C2'2  (bI  '3  '  b2'3'  b  1  'j  ')  H  C  3' 2  (  b  1 '2  '  b 2  '3  '  bl  '3  ')  f 


(36c) 


D3  =  r  jv3[i-  (V*')2]  T{Va  (b.vb2v-br2')  *c3'3(b,'2'b2'3'-b, v)[ 


(36d) 


DI  A'  j  C  1  h  (bt  '3  '  b2  V  bl'2')  ~‘Vl  1  “  (biv)2]  +<Vl  ^1  Vbl  V~b2'3')| 


(36e) 


D2  "A-  jCl'a(bi'3'b2'a'  S  v)  +C2'2  []  ~  (hiv)']  +  ^V2  (b.  VblV~  V3  )j 

(36f) 

D3  =  T'  |Cl'3  (bl  '3'b2'3'  -  b,  '2')  +  C2  3  [l~  (bi  v)2]  +C3'3  (bl'2'br3^-b2  '3')[ 


3'  1  \ 

D'  =  T'  )Cl'l  (bC2'b2'3'-bl'S')  +  C2'l  (bl'2'bl 


(36g) 


(36h) 
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°2  '  A'  |  C1  '2  (  b  I  '2  '  b2  '3  '  Vs')  +  C2  2  (blVblV  b2  '3  ')  +  c  3  '2  [ 1  "  (  b  1 '2  '  )  2  ]  | 


(36i) 


D3  ~  jCl'3  (blVb2V  Vs')  +  C2'3  (blVb|V  b2  '3  ')  +  C3  '3  [}  (  b  l  '2  ')  2  | 


L'  ~  l~  (bl'2')2  (b2'3')2  (blV  )2  +  2  blVbsVb 


1  3 


(36j) 

(36k) 


4.2  ELEMENTS  OF  THE  METRIC  TENSOR 

Expansion  of  Eq.  (35)  yields  the  expressions  for  the  rectangular  coordinates  as  a 
function  of  the  oblique  coordinates.  Differentiation  of  these  expressions  and  substitution 
into  Eq.  (29)  yield  values  for  the  elements  of  the  covariant  metric  tensor: 


gl  'l ' 

[(•..-) 

(••.•)*] 

(37a) 

S  2  'l  '  ~ 

[(«„•) 

(■>■•)  '] 

(37b) 

B»  V  - 

[(■»•) 

2*  ('2J')2* 

(*..■)*] 

(37c) 

*sV  = 

6  [  '2  '  = 

[M  (Cl20 

1  +  (car) 

(C22 ')  +  (C3 1  ') 

(C32')] 

(37d) 

63  V  = 

er3 '  = 

.(‘..'I  (Cl30 

+  (*=*»') 

(c23')  +  (cai') 

(C3  3  '  )] 

(37e) 

e3v  = 

&2  '3  '  = 

[(CL2')  (c13') 

1  -  ('»•) 

^C23 T  (C32') 

(c33 

(37f)) 

5.0  VECTOR  MAGNITUDE 


Expanding_^Eq.  (14  )  and  taking  the  square  root  yield  the  expression  for  the 
magnitude  of  V  in  terms  of  the  elements  of  the  covariant  metric  tensor  and  the 
contravariant  measure  numbers: 
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v  =  [e,v(",')!*eJv(v,')!-s,v(»s')2-%,v'1'v2' 

*2«i  vv'  v3  +  2  e2  vv!  v>  ]  (38) 

Solution  of  Eq.  (17)  for  the  v1  and  substitution  into  Eq.  (38)  give  the  expression  for 
vector  magnitude  in  terms  of  the  vector  components: 

V  =  [<V)*,+  (V)“,+  (V)J,+  2  g.v^iVBsv)  “  *  (V)t  -  (V)2  - 

+  2 Si '3 '(err 83  3')  ('V,VV  +  2  B2v(b2  V@3'3')  -'J  (VV(VV]  2  (39) 

The  general  expression  for  vector  magnitude  in  terms  of  vector  projections  is  quite 

complicated.  For  practical  applications,  however,  a  straightforward  procedure  is  to  solve 
.1  .* 

Eq.  (18)  for  the  vJ  in  terms  of  the  (VjJV  and  to  then  substitute  the  vJ  into  Eq.  (38). 

.1 

The  vJ  are  given  by 

v’  =  £  |(sih')lj  (Vi)  i'[b2'2'8»'3'-  (b2v)  2] 

-  ^82  '2  ')  2  (  ^  1  )  2  *  1  '2  '  £3  '3  '  —  81  '3  '  82  '3  '  J 

*(«.v)''('i)J{!r.-i,v-!,vt!y]  |  (40a) 

vS  "  5  j “  ( B 1  '> ')  '  (V1 )  1 '  [ei  '2 '  Sa's'  -  B'  s '] 

+  (e2v)'2  (Vi)z'  [srr83v-(siv)2  ] 

-  (8°  ->  •)  (Vl)  s'  [  8  ■  v  82  V  -  8 .  V  8,  v]  |  (40b) 

V”'  '  A  }(*>'■')  ‘''(Vi),'[8l'2'ea's'  ~8'VS!'!] 

-  (g2v)  14  (Va)2'[biT'B2'3'  -  B1VB1  v] 

'  (s.v)‘4  (vi)>-[s.VBsv-  (bi'«)‘]{  (40c) 


(40b) 
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where 

^  ”  Si 'l '^§2  V  S3 '3  '  —  (@2  's')  J  _  ®1  '2  '  ^81  '2  '  S3  '3  '  ~  Si  '3  '  g2  '3  '  J 

+  Si  V  [si's'-  ei'a'Srs']  (40d) 

6.0  ANGLE  SUBTENDED  BY  TWO  INTERSECTING  VECTORS 

Dividing  Eq.  (14)  by  V2  yields 

v1  V* 

1  =  6.V  T  —  (41) 

so  that  v*  /V  is  a  measure  number  of  the  unit  vector  u^.  Assume  that  a  vector  W 
intersects  V  to  define  the  angle  0v,w  as  shown  in  Fig.  3.  The  unit  vector  uw  and  unit 


Figure  3.  Angle  0vw_subtended  by  the  intersection 
of  vectors  V  and  W. 

t 

vector  measure  number  W1  / W  are  associated  with  the  vector  W.  The  magnitude  of  the 
vector  (uv  -  ifw )  can  be  found  by  the  cosine  law  to  be 


This  same  vector  magnitude  is  given  by  [see  Eq.  (24)] 


(42) 


(43) 


16 


AEDC-TR-79-35 


Comparison  of  Eqs.  (42)  and  (43)  shows  that 


arccos 


(44) 


6.1  ANGLE  SUBTENDED  BY  V  AND  A  RECTANGULAR  COORDINATE  AXIS 

The  set  of  measure  numbers  of  a  vector  is  a  univalent  (rank  one)  tensor  (Ref.  13.)  A 
contravariant,  univalent  tensor  is  defined  as  a  quantity  that  transforms  as  (Ref.  13) 

i 

U  "  rV  "  (45) 

t 

Once  again  let  x'  represent  a  rectangular  Cartesian  coordinate  system,  and  let  x1 
represent  an  oblique  Cartesian  coordinate  system.  A  unit  vector  Uj  in  the  x1  coordinate 
direction  has  a  single  non-zero  measure  number  ul  -  1.  The  corresponding  measure 
numbers  in  the  oblique  coordinate  system  are  given  by 


u 


cb1'  .  ( 

-  u  I  u 

J.  1  ' 


1  -  i,  UJ  =  U*  = 


0 ) 


(46) 


Now  Eq.  (44)  can  be  used  to  find  the  angle  subtended  by  the  vector  V  and  the  x> 
coordinate  axis. 


arc  cos 


(47) 


At  this  point  Eqs.  (36)  and  (45)  can  be  used  to  find  the 
u*  values  can  then  be  substituted  into  Eq.  (47  )  to  obtain 


u1 


in  the  terms  of  u1.  The 


^V.x* 


arccos 


(l  r  n1 '  l'  n2  2'  n3  3' 

}V  |elT'D,  V  +^‘2'2'D,  v  +  83VDj  v 
+  g1-„'Dr  vl  +  go'.-l)1  v2  i  e.’i'D1  v1 


S3 


2  -K  yl  +  8 2 ' 1  15 1  v"  '  8  ] 

’,2'+s5-2-d;  'v3']| 


(48) 


The  expressions  for  the  angle  0  in  terms  of  vector  components  can  be  found  by 
substituting  Eq.  (17)  into  Eq.  (48): 
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The  general,  3-space  expressions  for  <f>  in  terms  of  the  (V^),*  projections  are  quite 
unwieldy.  For  practical  applications,  the  measured  values  of  the  (VAV  can  be  inserted 
into  Eqs.  (40  a  through  c);  then  the  resultant  v‘  numerical  values  can  be  substituted  into 
Eq.  (48)  to  find  <p. 

7.0  COMPONENT  AND  PROJECTION  INTERRELATIONSHIPS 


For  purposes  of  analysis  it  is  often  convenient  to  use  expressions  that  interrelate 

,1 

components  and  projections  in  the  x1  and  xl  coordinate  systems.  The  familiar 
parallelogram  law  for  vector  addition  of  components  can  be  used  to  describe  the  vector 
quantity  V  in  the  x'  as  well  as  the  x‘  coordinate  systems  (see  Fig,  2): 


3 

V  =  L  (V)^ 

i=l 


(50) 
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and  3  ■ 

V=£(VV?,,  (51) 

l  =  t 

where  ui  and  Uj*  are  unit  vectors  collinear  with  the  x*  and  x*  axes  respectively. 

By  forming  the  scalar  product  ofV  as  expressed  in  Eq.  (50)  with  iq,  one  finds  that 


(V'ih 

1  0  o 

w 

on  2 

= 

0  10 

<v)2 

on3 

0  0  1 

-(V)3_ 

(52) 


which  simply  formalizes  the  earlier  statement  that  (Vj^)j  and  (V),  are  identical  quantities 
in  the  rectangular  Cartesian  coordinate  system  x'.  Forming  the  scalar  product  of  V  as 
expressed  in  Eq.  (51)  with  Uj'  yields  the  relationships  for  (VjJp  in  terms  of  (V)j-: 


O'j)  |  - 

1  b]V  b  ,  '3  ' 

(V)r 

OiV 

= 

b2'l'  1  b2  '3 ' 

cv)2 ' 

1 _ 

bn '  b3'2'  1 

^  O I3 

The  x1  components  in  terms  of  xl>  components  can  be  found  by 
product  of  Eq.  (50)  and  Eq.  (51)  with  u^  and  equating  the  results: 


(53) 

forming  the  scalar 


m,- 

r  cn-  c12'  c13- 

'  (v}l; 

O')  2 

ii 

c  21 '  c  22 '  c  23  ' 

(V)2- 

_<v 

c 3 1 '  c32 "  c33 '  . 

0  )3 

(54) 


The  functional  interrelationship  of  (Vi);-  and  (V),-  can  be  found  by  forming  the 
scalar  product  of  Eq.  (50)  and  Eq.  (51)  with  Uj1  and  equating  the  results. 
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t- 

- 

cn  ci'2  cr3 

1 - 

> 

(Vi)2- 

- 

c2'l  c2'2  c2'3 

(V)2 

.(VlV  . 

C3T  c3'2  c3'3 

(V')3  _ 

The  inverse  relationships  for  Eqs.  (52  through  55)  can  readily  be  found.  Because  of  its 
usefulness,  however,  the  inverse  relationship  for  Eq.  (55)  is  given  here: 


< 

^c2'2  c3'3  _  c3'2  c2'3^  ^c3 '2  c  1'3  ~  c  1'2  c3 '3^  1'2  c  2  '3  -  c  2'2  c  1'3^ 

’  (Vi)r' 

(V)2 

1 

=  A  • 

(c3T  c2'3  “  c2'l  c3'3^  I  '1  c3'3  _  c3'l  C1  '33  ^c2  T  c  1'3  ”  c  1  T  c 2 '3^ 

(Vl>2' 

CO 

> 

_  i 

^c2'l  c3'2  “  c3'l  c2'2^  ^c3'l  cl'2  ~  cl'l  c3'2^  l'l  c  2  '2  -  c2  'l  c  1  '2^ 

(Vi)3'_ 

(56) 


where 


A  = 


11  <c2'2  c3'3  c2'3  e  3  '2^  +  c  l'2^c  3  'l  c2'3  “  c2'l  c3'3^  +  cl'3(c2T  c3'2  ”  c3'l  c2'2^ 


8.0  2-SPACE  RELATIONSHIPS 


The  2-space  versions  of  the  preceding  relationships  are  important  because  of  the 
widespread  use  of  two-axis  LDV  systems  for  two-dimensional  flow  studies.  For 
simplification  let  x3  be  collinear  with  x2  ,  and  let  both  be  perpendicular  to  x1  ,x2  ,xl 
and  x2.  Then,  for  2-space  we  have  ■ 

v3  =  (V)3.  =  (^1)3 '  =  c3 1 '  =  c32 '  =  c3'l  =  c3'2 


“  D3  =  D3  =  D1  =  D2  =  D3  =  6lV  =  62V  =  0 


It  can  also  be  readily  shown  that 

c3'3  =  c33 '  =  8lV"  «2'2'“  1 
gl'2'  =  bl'2' 


D 


ln2y  [cri-c2Tbr2] 


(57) 

(58) 

(59) 

(60) 
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l2  =  sin  Cl'2  C2'2bl 


1  sin 


D2'  =  — L- 

2  siny 


n V  [“  '!'■  =2'.] 


[  ‘i-!lir‘'n] 


where  7  =  arccos  bjV- 


8.1  VECTOR  MAGNITUDE 


In  2-space,  Eq.  (38)  becomes 


V  ,  [(,■•)«.  (.2  )  2  -  2  b,  vJ  ",2”]  ' 


while  the  expression  for  V  in  terms  of  vector  components,  Eq.  (39),  reduces  to 


V  =  [(V)2,.(V)2.,2bi.2.(v)r(V)2,]l'2  (66) 

Substitution  of  Eqs.  (58)  and  (59)  into  Eqs.  (40a),  (40b),  and  (40d)  yields 

V’  =  [(Vl)|-kl’S'(Vi)s']  (67) 

v2  '  CA  [-br2-(vi)i-+  (v±)  2']  (68) 

Now  Eqs.  (67)  and  (68)  can  be  substituted  into  Eq.  (65)  to  find  the  simplified  2-space 
expression  for  V  in  terms  of  vector  projections: 

V  ■  jdvtfilK"  (Vf)  l-  2  b.v(vi),-(\)  =']['■  (69) 

It  can  be  seen  that  Eqs.  (66)  and  (69)  both  reduce  to  the  familiar  form  for 
rectangular  Cartesian  coordinates  when  7  =  90  deg' 

v  -  j  (v)  -  (V  \r  j1;  (70) 
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8.2  THE  ANGLE  <t> 


In  2-space,  the  simplified  expressions  for  0  in  terms  of  vector  components  are  found 


to  be 


1  =  arccos 


<0y  x2  =  arccos 


{■;[ci'i (VV]} 

(71) 

(72) 

The  expressions  for  0  in  terms  of  velocity  projections  are  given  by 

<0yxl  =  arccos]—!—  [  (c  ri  -  c2-j  bj  ^  )  (V  )  ^  +  (<=2'1  -  'll  bl  '2')(Vl)2  '1  [ 
|  Vsin  y  L  J  ) 


<0Vx2  =  arccos  —!—  (c  r2  ~  c2'2  b!  V  )  (Vx)  V  +  (c2'2  ”  C1'2  bl'2')  C V  j.)  2  '  |  [ 
I  VsirTy  -J  J 


It  can  be  seen  that  Eqs.  (71)  and  (73)  and  Eqs.  (72)  and  (74)  reduce  to  the  familiar 
forms  for  rectangular  Cartesian  coordinates  when  7  =  90  deg : 

■)  / 

<0V  xl  =  arccos  j  -  j 'j  CV)l'  +  c 2  'l  ( V ) 2  ' 


4>\j  X1  =  arccos 


|  y  [C1'2  (V)l  '  +  c2'2  (V^2' 


8.3  COMPONENTS  AND  PROJECTIONS 

In  2-space,  Eqs.  (52)  through  (56)  reduce  to 


1  0 

0  1 
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(78) 


(79) 


(80) 


(81) 


9.0  CONCLUSIONS 

The  complexity  of  the  general  3-space  relationships  clearly  illustrates  the  desirability 
of  implementing  orthogonal  LDV  measurement  axes  collinear  with  the  reference  system 
axes.  When  these  conditions  are  unattainable,  the  equations  in  this  report  can  be  used  to 
compute  the  values  for  velocity  magnitude  and  spatial  orientation  within  the  uncertainties 
dictated  by  the  remaining  LDV  system  error  mechanisms. 

The  provision  of  the  correct  equations  to  account  for  axis  nonorthogonality  and 
rotation  now  focuses  attention  upon  the  status,  with  regard  to  accuracy  and  precision,  of 
techniques  and  devices  for  determining  the  b,;j<  and  q'j  values  to  be  used  in  these 
equations.  Moreover,  further  reductions  in  LDV  instrument  errors  will  require  additional 
refinement  of  (1)  techniques  and  devices  for  determining  the  signal  frequency-to-veloeity 
calibration  factor  for  each  LDV  measurement  axis  and  (2)  techniques  and  devices  for 
measuring  the  signal  frequency  associated  with  each  axis. 

REFERENCES 

1.  Rabinawitz,  P.,  Jacobs,  S.,  Targ,  R.,  and  Gould,  G.  “Homodyne  Detection  of 
Phase-Modulated  Light.”  Proceedings  of  the  IRE,  Vol.  50,  No.  11,  November  1962, 
p.  2365. 


23 


AEDC-TR  79-35 


2.  Cummins,  H.  Z.  and  Knable,  N.  “Single  Sideband  Modulation  of  Coherent  Light  by 
Bragg  Reflection  from  Acoustical  Waves.”  Proceedings  of  the  IEEE,  Vol.  51,  No.  9, 
September  1963,  p.  1246. 

3.  Cummins,  H.Z.,  Knable,  N.,  and  Yeh,  Y.  “Observation  of  Diffusion  Broadening  of 
Rayleigh  Scattered  Light.”  Physical  Review  Letters,  Vol.  12,  No.  6,  10  February 

1964,  pp.  150-1 53. 

4.  Yeh,  Y.  and  Cummins,  H.  Z.  “Localized  Fluid  Flow  Measurements  with  an  He-Ne 
Laser  Spectrometer.”  Applied  Physics  Letters,  Vol.  4,  No.  10,  15  May  1964,  pp. 
176-178. 

5.  The  Accuracy  of  Flow  Measurements  by  Laser  Doppler  Methods.  Proceedings  of  the 
LDA-Symposium.  Copenhagen,  Denmark,  August  25-28,  1975. 

6.  Crosswy,  F.  L.  and  Brayton,  D.  B.  “Design  and  Evaluation  of  a  Laser  Doppler 
Velocimeter  for  the  AEDC-PWT  4-ft  Transonic  Tunnel.”  AEDC-TR-78-30,  July 
1978. 

7.  Rolfe,  E.,  Silk,  J.  K.,  Booth,  S.,  Meister,  K..,  and  Young,  R.M.  “Laser  Doppler 
Velocity  Instrument.”  NASA  CR-1199,  December  1968. 

8.  Barnett,  D.  O.,  and  Giel,  T.  V,  “Application  of  a  Two-Component  Bragg-Diffracted 
Laser  Velocimeter  to  Turbulence  Measurements  in  a  Subsonic  Jet.”  AEDC-TR-76-36 
(AD-A025355),  May  1976. 

9.  Stratton,  Julius  Adams.  Electromagnetic  Theory,  McGraw-Hill  Book  Company,  Inc., 
New  York,  194] . 

10.  Sokolnikoff,  Ivan  Stephen.  Tensor  Analysis.  Theory  and  Applications  to  Geometry 
and  Mechanics  of  Continua.  John  Wiley  and  Sons,  Inc.,  New  York,  1964. 

11.  Lass,  Harry.  Vector  and  Tensor  Analysis.  McGraw-Hill,  New  York,  1950. 

12.  Lennert,  A.  E.,  Brayton,  D.  B.,  and  Crosswy,  F.  L.  “Summary  Report  of  the 
Development  of  a  Laser  Velocimeter  to  be  Used  in  AEDC  Wind  Tunnels.” 
AEDC-TR-70-101  (AD871321),  July  1970. 

13.  Moon,  Parry  and  Spencer,  Domina  Eberle.  Vectors.  Van  Nostrand,  Princeton,  N.  J., 

1965. 


24 


AEDC-TR  79  36 


NOMENCLATURE 

Matrix  array  of  constants,  Eq.  (30) 
ap  Unitary  vector,  Eq.  (4) 

—y}  ^ 

a’  Reciprocal  unitary  vector,  Eq.  (6) 

bjy  Matrix  array  of  direction  cosines,  Eq.  (34) 

bi'1'  Inverse  of  bjy,  Eq.  (36a) 

cpj  Matrix  array  of  direction  cosines.  Eq.  (36a) 

Dj'  Matrix  array  of  constants,  Eq.  (36a) 

dr  Differential  displacement  vector,  Eq.  (3) 

ds  Differential  displacement  magnitude  in  the  x1  coordinate  system,  Eq.  (22) 

ds’  Differential  displacement  magnitude  in  the  x*  coordinate  system,  Eq.  (12) 

fj>  LDV  signal  frequency,  Eq.  (23) 

gjV  Covariant  metric  tensor,  Eq.  (12) 

gi V  Contravariant  metric  tensor.  Eq.  (13) 

ks  Scatter  radiation  propagation  vector,  Eq.  (23) 

ki  Illuminating  beam  propagation  vector,  Eqs.  (23)  and  (25) 

— ^ 

r  Position  vector,  Eq.  (3) 

u|  Unit  vector  collinear  with  the  x1  axis,  Eq.  (50) 

Iq'  Unit  vector  collinear  with  the  x1’  axis,  Eq.  (24) 

- — ^ 

O;  Unit  vector  collinear  with  V,  Eq.  (42) 

— T 

iitv  Unit  vector  collinear  with  W,  Eq.  (42) 

0  Volume,  Eq.  (5) 

V  Vector,  Eq.  (10) 
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Magnitude  of  V,  Eq.  (14) 

Component  of  V  collinear  with  the  x1  axis,  Eq.  (52) 

— y 

Orthogonal  projection  of  V  collinear  with  the  xL  axis,  Eq.  (52) 
Component  of  V  collinear  with  the  x1'  axis,  Eq.  (17) 

— y 

Orthogonal  projection  of  V  collinear  with  the  x»'  axis,  Eq.  (18) 

— ^ 

Covariant  measure  number  of  the  vector  V,  Eq.  (11) 

— y 

Contravariant  measure  number  of  the  vector  V,  Eq.  (10) 

Vector,  Eq.  (42) 

Magnitude  of  W,  Eq.  (43) 

Rectangular  Cartesian  coordinates,  Eq.  (2) 

Oblique  Cartesian  coordinates,  Eq.  (1) 

Oblique  Cartesian  coordinates.  Eq.  (31) 

Angle  equal  to  arccos  bj'2',  Eqs.  (60)  through  (64) 

A  constant,  Eq.  (56) 

A  constant,  Eqs.  (36b)  through  (36k) 

Kronecker  delta,  Eqs.  (7)  and  (19) 

Angle  subtended  by  LDV  transmitter  beams,  Eq.  (24) 

Wavelength  of  laser  radiation,  Eq.  (24) 

— y  — y 

Angle  subtended  by  the  intersection  of  V  and  W,  Eq.  (42) 

Angle  subtended  by  the  intersection  of  V  and  x*  coordinate  axis,  Eq.  (48) 
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